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Periodic multilayers of various periods were prepared according to an algorithm 
proposed by the authors. The reflectivity properties of these systems were inves- 
tigated using neutron reflectometry.The obtained experimental results were com- 
pared with the theoretical expectations. In first approximation, the results proved 
the main features of the theoretical predictions. These promising results initiate 
further research of such systems. 

1 Introduction 
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Neutron supermirrors are nowadays used in many physical experiments. They are mul- 
tilayer systems usually composed as a set of bilayers every one of which consists of two 
materials with high and low optical potentials. The supermirrors increase the angular 
or wave length range of total reflection comparing to mirrors consisting of a single ma- 
terial with high optical potential. If the single material gives total reflection for normal 
component k of the incident neutrons in the range < < /Cc, where kc is the limiting 
wave number for the given material, the supermirror can increase the interval up to nkc- 
A multilayer system that gives reflectivity 1 in the interval < A; < nkc is called Mn 
mirror. It became a usual practice to fabricate M2 and M3 mirrors. However there are 
also attempts to produce mirrors with higher n. The last record belongs to Japanese [1] 
who prepared mirror M6.7. Reflectivity of this mirror is shown in flg. [T} 

Last time all the mirrors were prepared in aperiodical fashion, which means that 
thicknesses of layers in bilayers vary with the bilayer number. F. Mezei and P. A. Dagleish 
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Figure 1: The measured neutron reflectivities of the Ni-Ti supermirrors fabricated with 
the large-scale ion-beam sputtering instrument (IBS) for the Spallation Neutron Source 
of J-PARC [1]. Reflectivities of 1) M3 mirror with 403 layers; 2) M4 mirror with 1201 
layers; 3) M6.7 mirror with 8001 layers. 
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performed the first experimental study of such supermirror in 1977 [2]. The algorithm for 
thickness variation was proposed by J.B. Hayter and H.A. Mook in 1989 [3]. According 
to this algorithm the change of thicknesses of neighboring bilayers is very small and does 
not match interatomic distance. It leads to creation of unavoidable roughnesses on layers 
interface. There exists also another algorithm proposed in [4], according to which the 
supermirror is to be produced as a set of periodic chains with some number N of identical 
bilayers. The variation of thicknesses of neighboring chains in this case is larger comparing 
to aperiodical systems, which may help to improve the quality of interfaces and therefore 
of the whole supermirror. 

The goal of the given work is to investigate how well we can control the thickness and 
quality of periodic multilayer systems prepared by Mirrotron Ltd, Budapest. In other 
words we want to see how well the neutron refiectivity of produced systems match theo- 
retical expectations, how large is diflFuse scattering because of technological imperfectness 
and whether we can explain and control them. 

Below we first present theoretical description of periodic multilayer systems, and calcu- 
late neutron refiectivities of periodic chains with N bilayers (N=2, 4, 8). Then we present 
results of measurements of refiectivities of fabricated multilayers and compare them to 
theoretical ones. 



2 Theoretical description of periodical chains 

A periodical chain consists of N bilayers every one of which contains two layers of diflFerent 
materials. In our case they were Ni and Ti. We call Ni with higher optical potential a 
"barrier", and accept the real part of its potential u^q^ which is 2.45 • 10"'' eV, for unity. 
It means that the other energies are defined in units of u^q. So, the full Ni potential with 
imaginary part is u^, = 1 — O.OOOMi. We call Ti with lower potential a "well", and its 
optical potential in units of ui^o is = —0.203 — 0.00012i. 

It was decided to investigate periodical stacks, that give Bragg refiection at the point 
k = 2. This point is the normal to the sample surface component of the incident neutrons 
wave vector, and its value is given in units of the critical wave number kc = ^/u^ of Ni. 
The point k = 2 have to be at the center of the Darwin table of the Bragg peak. Our 
main task is to find thicknesses of the two sublayers of a bilayer, to get Bragg peak (at 
oo) with maximal width of the Darwin table. 

Refiection amplitude of a periodical potential with N symmetrical periods is given by 
the equation [4] 

1 - exp{2iqNa) 
1 — it^ exp(2zgA/ a) 

where 

^ _ V(l + r)^-t^- V(l-r)^-t^ 
" V(l + r)2-t2 + ^(l_^)2_^2' 

V(l + t)2-r2 + ^(l-t)2-r2' 

and r, t are reflection and transmission amplitudes of a single period. 

In the case of a bilayer the potential of a period is not symmetric, as is shown in 
Fig. [2} Therefore we have to take into account a direction of reflection and transmission. 
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We denote 7^ the reflection amplitude for the incident wave propagating to the right, and 
V for the incident wave propagating to the left. Then 

-^ = n + tl-^^^ r = r^ + tl-^^^ t = (4) 
1 - rtr^ 1 - rtr^ 1 - r^r^ 



where 
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1 - ^06,^ exp(2zfc5,^/5,^) 1 - r^^^^ exp(2zfc5,^/5,^) 

rob,w ^ ^^^"^ - ut^^. (6) 

We see that the transmission amplitude, t, is symmetric, i.e. it is the same in both 
directions. We can also introduce the symmetrized reflection amplitude 

r = VT^V, (7) 

then, with account of asymmetry the equation ([T]) takes the form 

~R (h\-~R 1 - exp{2iqNa) 
^^^^^ - ^-l-RleM^^QNay 

where i?oo and Rn inherit the asymmetry of r, i.e. 

Rqo = —Rool Rn = —Rn^ (9) 

and i?Ar, Rqo are given by ([T]), ([2]) with symmetrized amplitude ([7| used for r. 

To flnd It^w of the layers in the bilayer, which at = 2 give the center of the widest 
possible Darwin table with |i?oo| = 1, we flrst neglect imaginary parts of the potentials 




Figure 2: A single element of a multilayered system is a bilayer composed of two differen 
materials. The layer of one material has a high optical potential u^, and a width Z^. The 
layer of the other material has lower potential and the width l^. 
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and represent t in the form t = \t \ exp^iip) with real phase ip. Then r = ±i\r \ exp{iijj) 
with the same phase ijj^ and Eq. ^ can be transformed to 

^ ^ VsinV;(fc) + |r(fc)| - ^smip{k) - \r{k)\ ^^^^ 
-^sin'0(/c) + \r{k)\ + ^sin'0(/c) — \r{k)\ 

From it we see that the Bragg reflection takes place when sin^ < |rp, the center of the 
Darwin table is at sin?/; = and the larger is |r|, the wider is the Darwin table. Therefore 
we must flnd the widths l^^^ from two conditions: 

iPih, k = 2) = 7r, |r(2)| = max(|r(/5, k = 2)|), (11) 

where we had shown dependence of ip and |r| on widths Z^^^. To have the conditional 
maximum at the point k = 2 we are to require maximum at this point of the function 

F{h, L) = \r{h, k = 2)\ + X[^Jj{h, A; = 2) - tt], (12) 

where A is the Lagrange multiplier. The maximum is found from three equations 

dli) dlyj dX ^ ^ ^'^ 



Solution of these equations gives 

kdu — k..,L., — 

2 

Therefore we must take 



kbh kJw = -• (14) 



= Z- = 0.907, = ^ JL = 0.766. (15) 

2^4 - Ub 2^4 - 

The unit of length, corresponds to fcc = 1 (Ni), therefore it is Ac/27r = 92 A. So = 83.4 
A, and = 70.4 A. It was decided to ask preparation of 3 samples with 2, 4 and 8 bilayers 
with thicknesses: Ni = 84 A, and Ti = 70 A. 




PSD 



Figure 3: The scheme of the neutron reflectometer at KFKI. The sample mirror can 
be rotated around an axis perpendicular to the plane of the Fig. The position sensitive 
detector os stationary and sufliciently large to register all the reflected neutrons, when 
grazing angle is sufliciently small according to experimental requirements. Collimation 
ande was 0.5 mrad and the neutron beam was monochromatic with wave length 4.28 A 
and presumable resolution ? 
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3 Measurement and processing of data 



In Figure [4] in linear scale are shown the results of fitting experimental data for 2, 4 and 
8 periods with the formula: 

R{k)= I \R^,(p)\^^ + n,, (16) 

k-d 

where 

Rns{p) = Rn{p) + TUp) ^ ^ (17) 

1 - RN{p)rso{p) 

is the reflection amplitude from a periodic chain of N periods evaporated over a semiinfinite 
substrate, and the substrate reflection amphtude is 

^ _ ^ 

rso = ^ _^ ks = \/k^ - Us. (18) 
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Figure 4: Fitting of the refiectivity data for periodic chains of bilayers evaporated on a 



thick fioat glass substrate. The fitting function is given by Equation (16). The data were 



obtained in Hungary at the refiectometer with wide stationary detector. The results are 
shown in linear scale for a) 2; b) 4; c)8 bilayers. The solid curves are theoretical curves 
with parameters found from fitting. 



Eq. (16) takes into account the final resolution of the installation, scheme of which 
is shown in Fig. [3| and possible existence of background, n^, in the system. The resolu- 
tion 6 and background were two fitting parameters. Besides them we took as fitting 
parameters the real parts of all the potentials u^^ u^^ the potential of substrate Ug in 
terms of that of Ni, and thicknesses of Ni and Ti layers in units of Ac/27r of Ni. Imaginary 
parts of the potentials were calculated from absorption cross sections to be: u^^ = 0.00014, 
u^^ = 0.00012, and u^^ = 0.0001. The last one was selected so small, because at first we 
thought that our substrate was pure silicon glass. The results for fitted parameters are 
presented in the first three lines of the Table |3| The last column of the table shows 
of the fitting. The pictures in Fig. [4] show a good fit of all the samples, however the 
thicknesses of the Ni layers are more than 20% higher and thicknesses of the Ti layers 
are more than 20% lower than in the project. The Ni potential was found to be slightly 
lower than is expected, which can be explained by presence of some oxygen or nitrogen 
impurities. The substrate potential was found to be too high comparing to pure silicon 
glass, but later we found that it was boron glass the potential Ug ^ OA for it is quite 
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0.964 


-0.258 


0.452 


0.00014 


0.00012 


0.0001 


1.121 


0.59 


0.036 


0.003 


25 


4 


0.934 


-0.388 


0.446 


0.00014 


0.00012 


0.0001 


1.182 


0.525 


0.033 


0.003 


114 
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0.993 


-0.242 


0.398 


0.00014 


0.00012 


0.0001 


1.061 


0.649 


0.035 


0.0089 


349 


8 


0.963 


-0.421 


0.415 


0.00014 


0.00012 


0.0001 


1.169 


0.54 


0.036 
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8 


0.972 


-0.349 


0.408 


0.00014 


0.00012 


0.0001 


1.13 


0.579 


0.036 


0.003 


151 



Table 1: Fitted values of real parts of potentials for Ni (b-barrier), Ti (w-well) and 
substrate (s), their thickness /, resolution 5^ background and for periodic chains with 
N=2,4,8 bilayers. Imaginary parts u^^ of the potentials were fixed. The results are given 
in dimensionless units. The unit of energy is equal to real part of Ni optical potential 
u^^- = 0.245 /xeV, and unit of length is the reduced critical wavelength Ac/27r = 92 A for 
Ni. The parameters at the 4-th line were obtained for zero background, and parameters 
of 5-th line were obtained for predetermined background = 0.003. It is interesting to 
see that the fitting with smaller number parameters (6 instead of 7) gave smaller x^- 




Figure 5: Logarithmic scale of diflFerent types of fitting of the refiectivity data from periodic 
chains of 8 bilayers on a thick fioat glass substrate, a) The fitting with 7 parameters, 
corresponding to Figure [i]^); b) Fitting with 6 parameters (n^ = 0); c) Fitting with 6 
parameters (n^ = 0.003); d) x-distribution in the case c). 

reasonable. The resolution 5 ^ 3.5 % and background = 0.003 are also quite accept- 
able. The worst was the value of the x^- It is especially high in the case of the 8 periods 
sample. 

The defects of fitting of the 8 periods sample are seen in Fig. [5] in logarithmic scale. In 
picture a), which corresponds to picture c) in Fig.jij it is seen that theoretical line goes too 
high at large k. It corresponds to overestimated background = 0.0089. If we exclude 
Til) from fitting parameters and put = the other fitting parameters are changed as 
shown in 4-th line of the Table |3j The result of fitting in logarithmic scale is shown in 
picture b) of the Fig. [5| It is seen that the background in this case is underestimated. If 
we put background = 0.003 fixed at the same level as was obtained for samples with 2 
and 4 periods, we obtain fitting parameters shown in 5-th line of the Table |3| and result 
of fitting in logarithmic scale shown in picture c) in Fig. [sj The parameter decreased 
(note that the number of fitting parameters in that case is 6, which is less than 7), however 
it is still too high, and in Fig. [sjl) there is presented the x distribution 

where y{kj) and (7{kj) are refiectivity and statistical error at experimentally measured 
points kj). This distribution has very high fiuctuations near minima of the data shown 
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Figure 6: Fitting of the glass reflectivity data, a) One fltting parameter u'] b) Two fltting 
parameters: u' and u"] c) Logarithmic scale of the Figure b); d) distribution for 

fltting with two parameters. 

at picture c) and near the potential edge. 

For analysis of the reason of so high fluctuations it was decided to analyze flrst the 
substrate. In Figure [6] there is presented the fltting of experimental data for substrate 
reflectivity fltted by function 

R{k) = |r,o(A;)P (20) 

with one fltting parameter u^^ while u^^ = 0.0002 was flxed. The result was = 0.405, 
which is in agreement with the fltting of periodic chains, and proves that the substrate 
is the boron glass, however = 170 was too high. In linear scale the graph is shown 
in Fig. [6^). The fltting with the same function but with 2 fltting parameters and u^^ 
gives = 0.406 and u^^ = 0.00464. In linear scale the graph is shown in Fig. ^) and 
in logarithmic scale in Figure [6]3). We see a strong deviation at large k. The parameter 

= 137 for such a fltting becomes a little bit less but remains still unacceptably high. 
The distribution x{k) over all the experimental points is shown in Figure |6]l). 

In Fig. g are shown the result of fltting with more parameters. The picture a) shows 
fltting with 4 parameters according to Equation 

k+S 

R{k)= I \ro{p)\'^ + n,. (21) 

k-S 

The fltting parameters were u^^ u" ^ 5 and u^j. It was obtained: u' = 0.408, u" = 0.003, 
S = 4% and = 0.0024, which is quite reasonable. The result of fltting in logarithmic 




Figure 7: Fitting of the data for glass reflectivity a) with 4 parameters and R{k) given 
by ( [21] ); b) Distribution of x{k) for such fltting. 



7 



scale is shown in Figure [T^). It looks quite well, however = 31, for this fitting is still 
too high. Distribution x{k)^ shown in Figure [7j3) demonstrates that there is still some 
peculiarity in experimental data near the critical edge. 



4 Scattering at the interface 

The anomaly near the potential edge of the substrate can appear because of scattering 
on surface roughness, on near surface inhomogeneities and even because of disordered 
distribution of atoms inside the glass [5]. The scattering at interface can be calculated 
with the help of distorted wave Born approximation (DWBA) i.e. with the help of Green 
function G{u^ fc, r, r^), which takes the interface into account. 



4.1 Green function of DWBA 

The Green function for refiection from the semiinfinite substrate satisfies the equation 

{A + P - e{z > 0)u) G{u, k, r, r') = -47r5(r - r^, (22) 

where u is the optical potential of an ideal medium at z > 0. Since the space along the 
interface is uniform, the Green function can be represented by two dimensional Fourier 
integral 



/d V 
-^exp {ip\\{r\\ - rj|)) G{u,p±, z, z'). 



Substitution of it into Eq. (22) shows that G{u,p±, z, z') satisfies the equation 
{(f/dz^ +pI- e{z > 0)u) G{u,p±, z, z') = 5{z - z'). 



(23) 



(24) 



According to common rules it can be constructed with the help of two linearly independent 
solutions '0i,2('^,P±, ^) of the homogeneous equation 



{dp /dz^ +p\- Q{z > 0)u) ^Iji^2{u,p±,z) = 0. 
For the function ipi{u^p±^ z) we can take 



(25) 



i)i{u,p^,z) = e{z < 0) 



exp{ik±z) + rso{k±) exp{—ik±z) 



+ Q{z > 0)tso{k±) exp{ik'j_z), 

(26) 

where rso{x) is refiection amphtude (18), and tsoix) = l+rso{x) is the refraction amphtude 



from vacuum into the substrate. For ^2(^1 P±, z) it is appropriate to take the function 

^2{u,p±, z) = e{z < 0) exp{-ip±z)to,{p±) + Q{z > 0) {exp{-ip'j^z) - r,o{p±) exp{ip'j^z)) , 

(27) 

where tos{x) = 1 — rso{x) is the refraction amphtude from substrate into the vacuum, and 
p'± = The function (27| contains incident wave propagating inside the matter. 

The Wronskian of the functions '0i,2(ii,P±, z) is Wi2(p±) = '2ip±tos{p±)- 



With functions ipi{u,p±, z) from (p6|) and ip2iu,p±, z) from (27) we have 



G{u,p±,z,z') = 



{Q{z > z')ipi{u,p±,z)'ilj2{u,p±,z') + Q{z' > z)ipi{u,p±,z')'ilj2{u,p±,z)) 



2ip±T'{p_i 



(28) 



4.2 Scattering because of disorder 

We first consider scattering because of disorder and incoherent scattering. Every atom j 
inside the ordered medium composed of coherent scatterers is enhghtened by the coherent 
wave field (f{rj) = tso{kj_) exp{ikrj) created by the incident wave exp(ifcr). However 
in presence of incoherency and disorder the enlightening field has fiuct nations, which we 
denote Sip^rj)^ and suppose that their average Sip{rj) = 0. These fiuctuations produce 
scattered field 

5^{r) = ^G{u,k,r,rj)5^{rj)bj = J d^p\\ex^{ipr)A{k,p), (29) 



3 

where 



A{k.p) = '-^^Y.^wHp'rj)Sip{rj)b,, (30) 

^ 3 

the wave vector p for the wave scattered into vacuum at z < is p = (pn, — and 



Pj^ = y/c^ ~ P^- suggest that scattering amplitudes bj of j-th atom, and its 

fiuctuating field are not correlated, therefore 8ip{rj)hj = 0, and the correlation of this 
product for different atoms is described by the correlation function 



S^{r,)b,5^{ri)k\ = Km^\\r,){S,i + g{\r, - n|)}, 

i.e. it is naturally supposed that correlation is proportional to |<yCi|^ itself. 

With the wave function (29) we can find the flux of scattered neutrons in the vacuum 
through any plane parallel to the substrate interface 

Ji = 



2i 



(5^*(r)-^(5^(r) - (5^(r)-^5**(r) 



(31) 

where * means complex conjugate, and S is some large area of the plane, over which we 
integrate. Ratio of this fiux to the incident one Sk± gives scattering probability 

wik) = ^ = ^ / f{\A{k,p)\^)d^p,. (32) 



Since for scattered waves p'^ = k'^^ and for propagating waves p± = y/c^ — is a real 
number, then 

d^P\\/p± = 2(PpS{p^ - k^) = 2^dp^, (33) 

and the last equality is correct when we can integrate over angle dcj) around normal. As 
a result (32) is transformed to 

k 

w{k) = ^ 1 ^|A(fe,p)|2dp^. (34) 



With it we can define differential probability, or indicatrix 

w{k^ ^ PA_) = i27rf^\A{k,pW, (35) 
9 



where 



\A{k,p)\^ 



tosiP±) 



27ipj 



(36) 



where k,^ = (fcy — py, /c^ + P±)' The double sum has diagonal part, which after transfor- 
mation to the integral over NqcPtj gives 



\A{k,p) 



2{k'i + pI) 



tso{P±)tsQ{k±) 



2ttp± 



and indicatrix 



Wd{k± p±) = No\b\^- 



2k±{k'i+p'[) ■ 

The nondiagonal part of the sum can be transformed to double integral: 
Y,F*irj)Fir,,) = No J d'r, J g{r,-r,,)F*{ri)F{r,,), 



(37) 



(38) 



(39) 



where 



g{r) = 5{r) - No-f{r) 



(40) 



is correlation function chosen in such a way as to exclude the diagonal part of the sum. 
For simphcity we shall take for 7(r) the Gaussian 



7(r) = (27r)-3/2 exp(-rV2a2), a = A^q"'^"- 



After substitution of (j){r) into (36) and integration we obtain the expression 

1 



K± 



l-exp{-a^K,^/2) 
pI + k'i 



■\tso{p±)tsoik±)\ 



,l-eM-Ap'±-k'A'm 



(41) 



(42) 



k±' p'i + k'i 

where in the last equality we put c = \bc\'^No, and approximated ~ [p'j_ — fc^]^, which 
is valid for small p'j_ and k'j_ and small angle between k and p. 






Figure 8: a) Function (44) in dependence on /c for c = 0.0001, u = 1 — «0.0002, a = 0.2, 
5 = 5; b) Function (43) in linear scale for several different p; c) the same as b) in 
logarithmic scale. 
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The function (42) represented in the form 

^na{k,p) = - eM-a'\p' - kf/2)], (43) 

and its integral, 

= c j - exp(-a>' - *'1V2)], (44) 



for c = 10~^, = 0.2, 6 = 5 are shown in Figure |8| The curves are sensitive to parameter 
a and the upper integration hmit b. With increase of b the left maximum in Figure 
increases comparing to the right one, and the right peak shifts to higher k. 

4.3 Fitting of the experimental data for boron glass 



The function (44) has two maxima, so we can hope to get better fitting for glass refiectivity 
shown in Figure [7^). In fitting we accepted potential of Boron glass to he u = 0.408 — 
0.003i as was obtained above. We also included smoothing of the interface with Debye- 
Waller factor, i.e instead of r^{k) we used r^(/c) = r^{k) exi[){—2h'^kk^)^ where 2/i^ was a 
fitting parameter. 

Scattering on randomness and fiuctuations was represented by function 

Mk) = c I | l'.°WUf)l^ |2 - exp(-aV' - kfm. (45) 





The number 2 in brackets shows that we take the sum of two functions (42) and (38). In 
the function (45) there are two fitting parameters: factor c and a^/2. 




1 10 
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Figure 9: a) Result of fitting in logarithmic scale of refiectivity from boron glass with 
function (|45]). The glass potential u = 0.408 — O.OOSz was fixed. The fitting parameters 
were found to be: c = 6.323 • 10-^ aV2 = 0.102, S = 0.022, rit = 2.401 • 10"^ and 
/i^ = 0.048. All the values are quite reasonable, b) x distribution for such fitting. We 
see that some fiuctuations still remain near the edge, but their amplitude decreased much 
comparing to that of Figure ItI: 



11 



The specular reflectivity was deflned as Rs{k) = |r^(/c)p — t^i(fc), because scattering 
decreases specular reflectivity. For the total reflectivity we used the expression 

R{k)= J R^{k)^+n, + w^{k), (46) 
k-d 

which contains two more fltting parameters: 6 and n^. Thus in total we had 5 fltting 
parameters. The result of fltting is shown in Figure [9j For upper limit of the integral 



6 = 9 in (|45|) the was 8.5. Though it is too high, it decreased considerably from value 
31. The most important result is that fltting shows a dip in reflectivity near the edge. 
Fitting parameters were found to be: c = 8.3 • 10~^, which corresponds to randomness of 
positions of all the atoms in the glass; = 0.204, which means that correlation length 
is near 40 A; resolution 5 = 2.2%, which is better than in previous flttings; background 
Til) = 2.4 • 10~^ is nearly the same as before; and 2h^ = 0.048, i.e. h ^ 15A, which means 
that in preparation of the float glass some Sn atoms diflFused into the glass to the depth 
of 15 A. Results of fltting shows that the dip near the edge is quite well described by 
nonspecular reflectivity, and is not a result of some instrumental peculiarity. The dip in 
X near the potential edge means that measured quantity is larger than theoretical one. It 
happens because in theory we does not take into account the neutrons scattered into the 
glass. These scattered neutrons after going through the sample can be registered by the 
detector. Because of them (the scattering is maximal near the edge due to Yoneda effect) 
experimental result is higher than theoretical one. 

4.4 Scattering on roughnesses at the interface. 

The usual approach is to consider roughness as a Gaussian process. It means that all the 
rough surface is treated as a single inhomogeneity and the scattered wave function is put 
down as 

^,(fc,r) = J d^piiexp{ipr)A{k,p), (47) 

where 

A(fe,p) = -iVo6^^^^^%M^;(,'). (48) 

Here 

C(r;,) 

/(k) = y J dz' exp{iKr') = ^ ^ exp(iK||rj|)[exp(mxC(r||)) - l]e(C(rj|) > 0), 



(49) 

if roughness is a cavity (note that its scattering density is —Nob, and k' = k'j_ + p'±), 
where q'j_ = \/ (f'^— uq), and 



A{k,p) = Nob J d'^r\ exp(iK||rj|) j dz'e{C{r\) < 0)x 

C(rj|) 

X [exp{ik±z') + rso{k±) exp{-ik±z')][exp{ip±z') + rso{p±) exp(-ip_L^')]) (50) 
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if roughness is a bump abov^the average interface. The parameter C is a random variable 
with probabihty density distribution 

P(C) = ^exp(-CV2a^), (51) 

where a characterizes the average height of roughnesses. Averaging both of expressions 
over ( we obtain corrections to specular reflectivity amplitude. It will be a combination 
of error function $(gcr), where q = k± ± p± or q = k^j_ + p^j_. 

For averaging of \A{k^p)\'^^ which depends on random variables ( at two different 
points rj| we need the density distribution of the Gaussian process 

where X is a correlation function, which depends on distance r\\ between two points, 
where (1^2 are deflned. 

Small roughnesses Though calculations with these formulas can be done up to the 
end without principal difliculties, we shall not proceed this complicated way and simplify 
our task assuming that the height, a, of roughnesses is sufliciently small [6], i.e. ak± <C 1. 
For small grazing angles it means that a <C 100 A, which is quite practical. In that case 
we can accept A{k^p) in the form 

A{k,p) = -iVo& '^'^^^.^l^"^^^^^ / dMar\^exp{^Kr^) (53) 
27Tk± J II II 

for all positive and negative (. With this function we do not have corrections to Tq, 
because {() = 0. The scattered waves are determined by 

j dM||exp(^^v||r;||) J (^M,, exp(^Av||r'2||)C(r;||)C(r'2||). (54) 

Averaging of C(^i||)C(^2||) ^^^^ (^2) gives 

(C(r';||)C(r''2||))=^'W||-r''2||). (55) 

Therefore 

{\Aik,p)\') = ^1^^ \tos{p±)tosik^)f j rfV^|exp(^Av||r^|)i^(r^|), (56) 

and indicatrix of nonspecular scattering is 

w{k ^ p) = 27r^^^ \tos{p)tos{k±)f j exp(^K||r^|)K(r^|). (57) 

To finish calculations we need to define correlation function. It is natural to suppose that 

K{r\\) = exp(-r^/2/2), (58) 



^In our geometry, where medium is at z > the bump is below the interface. 
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ay: 



0.2 0,4 0.S O.g 1 1.2 1.4 1.6 1.3 2 2.2 2.4 




.4 QJ 0.l5 0.7 O.g 1 1.1 IJ 1.3 1.4 IJ l.d 1.7 l.S 1.5 1 



Figure 10: a) Experimental data points and theoretical calculation of supermirror M2.b) 
Distribution It is obtained from direct comparison of experimental and theoretical 

date without fitting. 



where / is correlation length, or average dimension of roughnesses along the interface. 
With this function we have 



Kf{k\\) = J exp(zA^||r 11 )exp(-r J/2/2) = 27tI^ exp{-l^ 1^1/2), 



(59) 



and substitution into (57) gives 



w(k ^p) = \NM^^^ 

k 



(60) 



4.5 Scattering and fitting of periodic chains 

There are a lot of opportunities how to include roughness at interfaces in multilayer 
systems. We can suppose that roughnesses are independent on every interface, or they 
can correlate between interfaces. It seems, that with sufficiently many fitting parameters 
it is possible to fit any result. However we shall not go this way. Because of so many 
opportunities it is better first to study experimentally the angular distribution of non 
specularly refiected neutrons, find its distinctive features and after that compare them 
with theoretical predictions based on different theoretical models. This is the way we are 
going proceed further. 



5 Supermirror 

Besides of the samples described above there was also prepared in BNC Budapest a 
supermirror M2, which consisted of 8 periodic chains and total number of 59 bilayers. 
The periods and number of them were found according to the above prescription with 
some corrections. The result of measurements comparing to calculations, in which we put 



Ni potential to be 0.964-0.0051 and Ti potential -0. 26-0.0051, is shown in Figure [10^). We 
see good coincidence. If we limit calculation of to the range /c < 2, then = 7, which 
is not bad, if to take into account that there were no fitting at all except some guess 
about imaginary parts of the potentials. The imaginary parts are higher than table ones 
because of possible impurities, inhomogeneities and surface roughness. 
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6 Conclusion 

Cooperation of theoreticians and experimentalists in research of multilayer systems is 
found to be very fruitful. We see that technology of preparation of such systems by 
Mirrotron Ltd, Budapest is good, but it can be further improved after analysis of surface 
imperfection and their correlation with parameters of producing systems. This analysis 
can be performed with new experiments aimed at investigation of diffuse scattering and 
angular distribution of reflected neutron with better angular resolution. 
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